Tunneling spectra of strongly coupled superconductors: Role of dimensionality 



C. Berthod 

DPMC-MaNEP, Universite de Geneve, 24 quai Ernest-Ansermet, 1211 Geneve 4, Switzerland 

(Dated: July 13, 2010) 

We investigate numerically the signatures of collective modes in the tunneling spectra of superconductors. 
The larger strength of the signatures observed in the high-r f superconductors, as compared to classical low-r c 
materials, is explained by the low dimensionality of these layered compounds. We also show that the strong- 
coupling structures are dips (zeros in the d 2 I/dV 2 spectrum) in rf-wave superconductors, rather than the steps 
(peaks in d~I /dV~) observed in classical s-wave superconductors. Finally we question the usefulness of effective 
density of states models for the analysis of tunneling data in <i-wave superconductors. 
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I. INTRODUCTION 

Many experiments have shown that the electrons in cuprate 
high-T ( superconductors (HTS) are significantly renormalized 
by the interaction with collective modes. This renormal- 
ization appears in photoemission measurements as velocity 
changes in the quasi-particle dispersion (the "kinks") accom- 
panied by a drop of the quasi-particle life-time. 1 ' 2 In tunnel- 
ing, the renormalization shows up as a depression, or "dip", 
in the dl/dV curve with the associated nearby accumula- 
tion of spectral weight (the "hump"). 3 Similar signatures ob- 
served by tunneling spectroscopy in classical superconductors 
were successfully explained by the strong-coupling theory of 
superconductivity. 4-6 There are, however, two striking differ- 
ences between the structures observed in the cuprates and in 
low-T c metals such as Pb or Hg. The dip in the cuprates is 
electron-hole asymmetric, being strongest at negative bias, 
while no such asymmetry is seen in lead. The electron-hole 
asymmetry of the dip is due to the electron-hole asymmetry 
of the underlying electronic density of states (DOS). 7 ' 8 Sec- 
ond, while the structures are subtle in low-r t materials — they 
induce a change smaller than 5% in the tunneling spectrum — 
the cuprate dip is generally a strong effect which, for instance, 
can reach 20% in optimally doped Bi2Sr2Ca2Cu30 lu+ ,5 (Bi- 
2223). 9 It is tempting to attribute this difference of intensi- 
ties to a difference in the overall coupling strength, as sug- 
gested by the largely different T c values. However, a compari- 
son of the effective masses indicates that the couplings are not 
very different in Pb where 10 m*/m = 2.1 and in the Bi-based 
cuprates where m* /m varies between 1.5 and 3 as a function 
of doping. 1112 Here we show that the large magnitude of the 
dip feature results from the low dimensionality of the materi- 
als and the associated singularities in the electronic DOS. 

Tunneling experiments in strongly coupled classical super- 
conductors have been interpreted using a formalism 13 that ne- 
glects the momentum dependence of the Eliashberg functions 
and of the tunneling matrix element, and further assumes that 
the normal-state DOS iVo ( (0 ) is constant over the energy range 
of interest. The tunneling conductance, then, only depends 
on the gap function A(a>), whose energy variation reflects the 
singularities of the phonon spectrum. 1415 The dimensional- 
ity of the materials does not enter in this formalism. The ef- 
fect of a non-constant Nq(co) on the gap function has been 
discussed in the context of the A15 compounds. 16 However 



the direct effect of a rapidly varying Nq((o) on the tunnel- 
ing conductance became apparent only recently in the high-r c 
compounds, 817 ' 18 and requires to go beyond the formalisms of 
Refs 13 and 16. In particular, one can no longer assume that 
the tunneling conductance is proportional to the product of the 
normal-state DOS by the "effective superconducting DOS" 13 
Re [|o|/-\/© 2 — A 2 (o)J , so that nothing justifies a priori to 
normalize the low-temperature tunneling conductance by the 
normal-state conductance as was done with low-T c supercon- 
ductors. 

Among the new approaches introduced to study strong- 
coupling effects in HTS, some have focused on generalizing 
the classical formalism to the case of of-wave pairing, 19-22 
still overlooking the dimensionality. Other models are strictly 
two dimensional (2D) and pay attention to the full electron 
dispersion, 7 ' 8 ' 17 ' 23-25 taking into account the singularities of 
Nq(co). Most of these studies assume that the collective mode 
responsible for the strong-coupling signatures is the sharp 
(tt, 7t) spin resonance common to all cuprates near 30-50 meV 
(Ref. 26), but a phonon scenario was also put forward. 22 ' 23 In 
the present work, we extend these approaches to three dimen- 
sions (3D) by means of an additional hopping f, describing 
the dispersion along the c axis, and we study the evolution of 
the strong-coupling features in the tunneling spectrum along 
the 2D to 3D transition on increasing t z . For simplicity we 
restrict to the spin-resonance scenario; the electron-phonon 
model can be treated along the same lines, and both mod- 
els lead to the same main conclusions. The model we use 
is described in Sec. II, results are presented and discussed in 
Sec. Ill, and Sec. IV is devoted to investigating the validity 
and usefulness of the effective superconducting DOS concept. 



II. MODEL AND METHOD 

Following previous works 817 ' 27 we assume that the differ- 
ential conductance measured by a scanning tunneling micro- 
scope (STM) is proportional to the thermally-broadened lo- 
cal density of states (LDOS) at the tip apex, 28 ' 29 and that fur- 
thermore the energy dependence of the LDOS just outside the 
sample follows the energy dependence of the bulk DOS N(co): 



^« / dco[-f'(co-eV)]N(co), 



(1) 



2 



where /' is the derivative of the Fermi function. In STM ex- 
periments, various sources of noise may contribute to broaden 
N(co) further; when comparing theory and experiment we 
shall take these into account by a phenomenological Gaussian 
broadening. 

In the superconducting state, the interaction with longitudi- 
nal spin fluctuations is described by the 2x2 Nambu matrix 
self-energy, 



t(k,(o) = -j- J £^g 2 Xs (q,i£ln)x 

Goik-qJCOn-iSln)].^ 



•co+iO 4 



(2) 



with Xs the (S Z S Z ) spin susceptibility, g = \/3hJ/2 the cou- 
pling strength with / the spin-spin interaction energy, and Go 
the 2x2 BCS matrix Green's function. The sums extend over 
the jY vectors q = (q^ ,q z ) in the three-dimensional Brillouin 
zone and the even Matsubara frequencies z'£2„ = 2nn/j3 with 
j3 = (IcbT)~ 1 . Equation (2) gives the lowest-order term of an 
expansion in 7. 30 We follow Ref. 7 and use for Xs a model in- 
spired by neutron-scattering experiments on the high-r c com- 
pounds. In this model, Xs has no q z dependence and is the 
product of a Lorentzian peak centered at (n, n) in the 2D Bril- 
louin zone and another Lorentzian peak centered at the reso- 
nance energy i2 sr . The widths of the peaks are Aq in momen- 
tum space and r sr in energy. This simple separable form of Xs 
allows to evaluate analytically the frequency sum in Eq. (2), 
and to perform analytically the continuation from the odd fre- 
quencies ico„ = (2n + to the real-frequency axis. The 
remaining momentum integral is a convolution which can be 
efficiently performed using fast Fourier transforms. Also, the 
absence of q z dependence in Xs implies that the self-energy 
does not depend on k z . We use a BCS Green's function broad- 
ened by a small phenomenological scattering rate F, 



G (k,ico„) 



[ico„ + lT(iOn)]% + &T3 + A k ti 



[i(D n + iT(i(o n )] 2 -&- 
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where F(ico„) = rsign(Im/a)„), f, are the Pauli matrices with 
% the identity, and Ajt = Ao(cosk x — cosk y ) /2 is the c/-wave 
gap which we assume k z independent for simplicity. The addi- 
tional effects resulting from a possible weak modulation 35 - 36 
of the gap along k z will be discussed toward the end of Sec. III. 
We do not address here the origin of the pairing leading to the 
BCS gap Ak . With the high-r c - compounds in mind, we con- 
sider a one-band model of quasi-2D electrons with a normal- 
state dispersion, 



& = + 2t z cos(k z c), 



(4) 



where §j = £j — jU, fl being the chemical potential and 
e^H a five-neighbor tight-binding model on the square lat- 
tice (fl = 1), & = 2fi (cOSfc v + COSfey) + 4t2COSk x CO&k y + 

2tj(cos2k x + cos2k y ) + 4?4(cos2A: v cosA: v + cosk x cos2k y ) + 
4?5 cos 2k x cos2k y . 

The momentum dependence of the self-energy in Eq. (2) is 
not small (see, e.g., Fig. 1 of Ref. 7 and Fig. 2 below). This is a 



major difference with respect to the electron-phonon models 
describing low-r c three-dimensional metals, where the mo- 
mentum dependence of the self-energy can be neglected. The 
calculation of the DOS is therefore much more demanding 
since two three-dimensional momentum integrations must be 
performed for every energy co. The DOS is given by 



N(co) 



Af H ( TtJ 



ImG n (fc,Ct)), 



(5) 



where G\\ is the first component of the matrix G(k,Co) = 
[GQ l (k,(o) — £(Jt,o))] _1 . In Eq. (5), the k z integration can 
be performed analytically (see Appendix) but not in Eq. (2). 
In order to achieve a good accuracy when computing the DOS 
N{co), we evaluate the self-energy using a 2048 x 2048 x 256 
mesh in momentum space and a value T = 0.5 meV. For 
the evaluation of the tunneling conductance T is increased 
to 2 meV, which allows to decrease the mesh size to 1024 x 
1024 x 256. 

The model in Eqs. (l)-(5) has several parameters but our 
focus is on the c-axis hopping energy t z . In HTS, t z is not 
larger than a few meV, and setting it to zero seems appropri- 
ate to discuss tunneling data. Indeed, in the 2D limit the model 
was found to fit the experimental data for optimally doped Bi- 
2223 very well. 8 27 Here we take the parameters determined 
from one such fit as a starting point, and we vary t z to demon- 
strate the role of the dimensionality on the tunneling spec- 
trum. The band parameters t\ 5 are —200, 31, —16, 8, and 
—7 meV, and the chemical potential is jj. = —200 meV. The 
gap magnitude is Ao = 46 meV. The spin-resonance energy 
is £2 sr = 34 meV, its energy width r sr = 2 meV, and its mo- 
mentum width Aq = 1.6a . Finally the coupling strength is 
g — 775 meV, which implies a quasi-particle residue Z = 0.44 
and a mass renormalization m* / m = 2.32 at the nodal point of 
the Fermi surface. The temperature is set to T = 2 K unless 
stated otherwise. The resulting theoretical tunneling conduc- 
tance for t z = is compared with experimental data in Fig. lc 
(topmost curve). 



III. RESULTS 

The evolution of N((o) with increasing t z is displayed in 
Fig. la. In the 2D limit, the DOS shows sharp and particle- 
hole asymmetric coherence peaks, strong and asymmetric 
dips, as well as humps and shoulders where the spectral 
weight expelled from the dips is accumulated. This produces, 
in particular, a characteristic widening of the coherence peaks 
basis, which becomes triangular. The particle-hole asymme- 
tries reflect the particle-hole asymmetry of the corresponding 
bare DOS No(co) shown in Fig. lb, whose Van Hove singu- 
larity (VHS) lies slightly below Ep at —16 meV: on the one 
hand, the spectral weight of the VHS goes to a larger de- 
gree into the negative-energy coherence peak, and on the other 
hand the enhancement of the scattering rate due to the VHS 
is stronger at negative energy, explaining the stronger dip at 
(0 < 0. 7 ' 8 This can also be seen in Fig. 2 where the electron- 
scattering rate — ImEn(fc, co) is displayed for the nodal and 
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FIG. 1. (a) Density of states N(co) for superconducting electrons cou- 
pled to a (n, %) spin resonance, as a function of the c-axis hopping 
t z . For r, = 0, the system is two dimensional while for t z = 200 meV, 
it is three dimensional, (b) Normal-state bare DOS Nq((0) for the 
same t z values, (c) Tunneling conductance for the same t, values. 
The temperature is T = 2 K and a Gaussian broadening of 4 meV 
has been applied. The circles show the experimental data of Ref. 8 
compared to the t- = spectrum. The curves in (a), (b), and (c) are 
offset vertically for clarity. 
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FIG. 2. Scattering rate — ImZn as a function of energy for several 
values of t z at (a) the nodal point k\\ « (0.41,0.41);r/a and (b) the 
and nodal point Jtii ~ (l,0.05)7r/a of the Fermi surface shown in the 
insets. Curves are offset vertically for clarity. The dashed vertical 
lines delimit the energy range |co| < fi sr where inelastic scattering 
by the spin resonance is forbidden. The dotted vertical lines indicate 
±(£i sr + A ). 



anti nodal points of the Fermi surface. The scattering rate van- 
ishes for \ co\ < £2 S1 - and has a pronounced, particle-hole asym- 
metric maximum near | (0 1 = n sr + Ao (more precisely between 
£2 sr + Ao and £2 sr + ( ^ + A§] 1//2 ). It is also clear from the 
figure that the energy of the scattering-rate peak shows no dis- 
persion with momentum 7 but its intensity is strongly momen- 
tum dependent and larger by a factor ^2.5 in the anti nodal 
region as compared to the nodal region. 

For t z > 0, the logarithmic divergence in No{co) is cut on 
the scale of 4t z due to dispersion along the c axis (Fig. lb). 
No significant change in either N(co) or dl/dV is observed 
for t z = 10 meV. This value is an upper bound for the c-axis 
hopping energy in the cuprates, and the relative insensitivity 
of the DOS to a small c-axis dispersion justifies the use of 
two-dimensional models for these systems. At larger t z val- 
ues, however, the suppression of the divergence in Nq(co) in- 
duces a drop of the coherence peaks in N(a>) and dl/dV. This 
is a direct effect of dimensionality on the tunneling spectrum, 
which was overlooked in the conventional strong-coupling ap- 
proaches of Refs 13 and 16. Simultaneously the peak in the 
scattering rate is also suppressed with increasing t z (Fig. 2), 
leading to a weakening of the dip feature in N(co) and dl/dV. 
This is an indirect effect of dimensionality, that is only re- 
vealed in the strong-coupling signatures. 

As Fig. 2 shows, increasing the dimension not only sup- 
presses the peak at i2 sr + Ao in the scattering rate but it also 
reduces its momentum dependence. In 2D, this peak arises 
because the q sum in Eq. (2) is dominated by the saddle- 
point region near &m = (71,0) and jfc M / = (0,7T), where the 
spectral weight of the BCS Green's function is largest — 
i.e., £|| — <7|| ~ t MM i. Hence the peak energy is determined 
chiefly by the BCS excitation energy at k M M /, shifted by i2 sr 



due to the convolution with the spin susceptibility, and the 
peak intensity is controlled by the momentum dependence of 
Xs (k — k MM /,Q. SI ), which is at maximum for kit = k M > M . 
The situation changes in 3D because the anti-nodal regions no 
longer dominate the spectral weight, as illustrated in Fig. 3. 
This figure displays the partial BCS density of states, i.e., the 
part of the BCS DOS originating from states close to the (n, 0) 
and equivalent points. While in the 2D limit, a region cover- 
ing just 14% of the zone around (71,0) provides 56% of the 
spectral weight for energies between i2 sr and £2 sr + Ao, its 
contribution is reduced to 21% in the 3D limit. Hence the 
scattering rate in 3D is nearly momentum independent and al- 
most constant above i2 sr + Aq. Finally, the 2D to 3D transition 




FIG. 3. Partial BCS density of states in two (2D) and three (3D) 
dimensions. The thin solid lines show the contribution coming from 
the (%, 0) region of the two-dimensional Brillouin zone, shaded in 
black in the inset while the dashed lines show the contribution of the 
remainder of the zone (shaded in gray). The thick line is the total 
BCS DOS. 
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also suppresses the particle-hole asymmetry of the scattering 
rate. This again results from the disappearance of particle- 
hole asymmetry in the underlying bare DOS (Fig. lb) and 
in the corresponding BCS DOS (Fig. 3). Thus the k z dis- 
persion simultaneously defeats four players who contribute to 
make the strong-coupling signatures in the 2D high-r c - super- 
conductors distinctly different from those in 3D metals: the 
Van Hove singularity, the particle-hole asymmetry, the mo- 
mentum dependence, and the strong scattering enhancement 
at |fi)| w £2 sr + Ao, especially near (n,0). 

In the curves of Figs, la and lc corresponding to the 3D 
limit, the strong-coupling signatures are barely visible. Their 
magnitude is ~ 1%, smaller than the ~ 5% value observed 
in Pb. The origin of this difference lies in the gap symmetry. 
In J-wave superconductors, the coherence peaks in the BCS 
DOS are weak logarithmic singularities 37 while in s-wave su- 
perconductors, they are strong square-root divergences. The 
strength of the scattering-rate peak at Q. SI + Ao, and conse- 
quently the strength of the dip in the DOS and tunneling 
spectrum, are determined by the strength of the coherence 
peaks in the BCS DOS, as is clear from Eq. (2). In the case 
of a cZ-wave superconductor, the coherence peaks are cut in 
3D as compared to 2D (see Fig. 3) in the same way as the 
logarithmic VHS in Fig. lb, resulting in the suppression of 
the scattering-rate enhancement at £2 sr + Ao in Fig. 2. (Note 
that, roughly speaking, the scattering rate is proportional to 
the BCS DOS shifted in energy by ±£2 sr .) The suppression 
of the BCS coherence peaks with increasing dimension also 
occurs in .s-wave superconductors, but with one difference: 
if, on the one hand, the part of the coherence-peak spectral 
weight coming from the VHS gets suppressed, on the other 
hand, the square-root gap-edge singularities persist in any di- 
mension. Therefore, in .s-wave superconductors the strong- 
coupling signatures remain clearly visible in 3D. This is il- 
lustrated in Fig. 4a. The 2D and 3D DOS curves of Fig. la 
are compared to the curves obtained for the corresponding s- 
wave model, i.e., with all parameters unchanged except the 
gap which is replaced by A^ = Ao = 46 meV The changes 
are quite dramatic. The first effect to notice is a drastic reduc- 
tion in the peak-to-peak gap A p in the .s-wave case: a conse- 
quence of the pair-breaking nature of the coupling Eq. (2) in 
the .s-wave channel. 38 ' 39 Still, the strong-coupling signatures 
appear at the same energy £2 sr + Ao = 80 meV in both d and s 
wave, due to our choice of the lowest-order model E « % s Gq 
in Eq. (2). The second observation is that the strong-coupling 
signatures look like steps in the .s-wave DOS, like in the clas- 
sical superconductors, 13 reflecting the asymmetric shape of 
the BCS .s-wave coherence peaks. In contrast, the signatures 
appear as local minima in the c/-wave DOS, because the co- 
herence peaks of the J-wave BCS DOS are nearly symmetric 
about their maximum. In short, the strong-coupling features 
give an "inverted image" of the BCS coherence peaks. 8 An in- 
teresting consequence follows: while in .s-wave superconduc- 
tors, the strong-coupling structures correspond to peaks in the 
second-derivative d 2 I/dV 2 spectrum, for a of-wave gap they 
correspond to zeros in the d 2 I/dV 2 spectrum, as demonstrated 
in Fig. 4b. This conclusion applies equally to phonon mod- 
els and calls for a reinterpretation of cuprate d 2 I/dV 2 data in 
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FIG. 4. (a) Comparison of the d-wave and .?-wave DOS for t z = 
(2D) and t z = 200 meV (3D). The thick lines show N(m) as in Fig. la. 
The thin lines show N(co) computed with an .s-wave gap of the same 
magnitude Ao = 46 meV, and all other parameters unchanged, (b) 
Second-derivative tunneling spectrum, d 2 I/dV 2 , in the region of the 
strong-coupling signature at negative bias. In the .s-wave case, there 
is a peak in d 2 I/dV 2 close to the energy £2 S1 - + Ao while in the rf-wave 
case, there is a sign change in 2D and no clear signature in 3D. 



which peaks were assigned to phonon modes. 40 ' 41 Finally, one 
sees from Fig. 4 that in 3D the signatures remain strong for an 
s-wave gap, for the reason explained above, while they have 
almost disappeared in the of-wave case. 

The previous discussion underlines the role of the BCS 
coherence peaks in the formation, strength, and shape of 
the strong-coupling signatures. More generally, for such 
signatures to occur there must be divergences (or at least 
pronounced maxima) in the non-interacting DOS. Peaks in 
the "bosonic" spectrum are not sufficient, although they are 
necessary. Indeed, phonon structures are absent from the 
normal-state spectra of classical superconductors 32 because 
the normal-state DOS is flat, in spite of the facts that the 
phonon spectrum and the electron-phonon coupling do not 
change significantly at T c . In contrast, the normal-state DOS 
of 2D high-r c superconductors exhibits structures, either the 
pseudogap 3 or the bare VHS. 18 One can therefore expect to 
see strong-coupling features in the normal-state spectra of 
HTS, provided that the peaks in the bosonic spectrum sub- 
sist above T c . Figure 5 (thin lines) shows the normal-state 
DOS implied by setting Ao = in our model, keeping the 
other parameters fixed (including temperature). As expected 
sharp strong-coupling features remain in 2D at energies ±i2 sr 
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-50 meV while nothing but very weak 



structures subsist in 3D, signaling the onset of scattering at 
±n sr . Unfortunately it turns out that in the HTS the spin reso- 
nance is absent above T c — or at least below the background 
level of neutron-scattering experiments. 42 The normal state 
of Bi-2223 has not been investigated by neutron scattering 
so far but we may borrow information from the much stud- 
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FIG. 5. Normal-state DOS. The thin lines show the T = 2 K DOS 
for A = 0; the thick lines show the T = 200 K DOS for A = and 
r sr = 14 meV (see text). All other parameters are as in Fig. la. 



ied YBa2Cu306 +x system (Y-123). In Y-123, the normal-state 
spin susceptibility preserves its separable form with indepen- 
dent momentum and energy variations. 43 It is still centered at 
(tt, n) with a broad maximum at a characteristic temperature- 
dependent frequency i2 s f w i2 sr . For the purpose of illus- 
trating the effect of a broad spin-fluctuations continuum on 
the normal-state tunneling spectrum, it is sufficient to use the 
same model as in the superconducting state but with the new 
parameter r sr = 14 meV. 44 The resulting DOS calculated at 
T = 200 K is shown by the thick lines in Fig. 5. The strong- 
coupling signatures are almost washed out in 2D and com- 
pletely in 3D. This is not due to the thermal broadening of 
Eq. (1), not included in the DOS N(co), but mostly to the in- 
trinsic temperature dependence of the self-energy in Eq. (2), 
and, to a lesser extent, to the broader spin response. Hence, if 
structures due to interaction with spin fluctuations are unlikely 
to show up in the normal state of HTS, those associated with 
the interaction with phonons may well be observable if the 
coupling is strong enough since this coupling will not change 
appreciably at T c . 

In the present study, we have overlooked a possible k z de- 
pendence of the BCS gap, retaining only the k z dependence 
of the bare dispersion. A weak modulation of the BCS gap 
along k z is expected in 3D systems. 35 As shown in Ref. 36, 
such a modulation has the effect of cutting the logarithmic co- 
herence peaks on the scale of 2A Z , with A z the amplitude of 
the gap modulation. This is similar to the effect of t z on the 
BCS coherence peaks, which are cut on a scale correspond- 
ing to the gap variation along the warped 3D Fermi surface, 
namely, ~ Aot z /4t\, as seen in Fig. 3. The expected effect of 
A- on the scattering rate is also an additional broadening on 
top of the one produced by t z , Gq being replaced by its k z av- 
erage in Eq. (2). Therefore, we expect that the gap modulation 
along k z will contribute to suppress the coherence peaks and 
the strong-coupling features even further with increasing t z , as 
compared to the results in Fig. 1. 

Our results can be summarized as follows. The formation 
of clear strong-coupling structures in the tunneling conduc- 
tance requires two ingredients: (A) at least one peak in the 
spectrum of collective excitations and (B) at least one peak 
in the non-interacting or superconducting DOS. In classical 
superconductors, (A) is provided by optical phonons and (B) 
is the asymmetric square-root singularity at the edge of the 
i-wave gap: strong-coupling features are asymmetric steps — 



peaks in the d 2 I/dV 2 curve — and dimensionality plays no big 
role because (A) and (B) are present in any dimension. In 
the normal state, there is no signature because (B) is absent. 
In high-r c layered superconductors, (A) is provided by the 
spin resonance and (B) has two sources: (Bl) the logarith- 
mic Van Hove singularity in the bare DOS; (B2) the symmet- 
ric logarithmic singularities at the edge of the of-wave gap. 
Strong-coupling signatures appear as local minima — zeros in 
the d 2 I/dV 2 curve — but they vanish with increasing dimen- 
sionality from 2D to 3D because (Bl) and (B2) both get sup- 
pressed by the c-axis dispersion. In the normal state of two- 
dimensional HTS, (B2) is absent, leaving aside the question 
of the pseudogap but (Bl) remains and strong-coupling sig- 
natures are thus expected unless (A) disappears at T c . This is 
the case for the spin resonance but certainly not for phonons, 
leaving open the possibility that phonon structures might be 
observable in the normal-state tunneling spectra. 



IV. DOS AND EFFECTIVE DOS 

The conventional theory of electron tunneling into 
superconductors 13 leads to an equation identical to Eq. (1) 
for the tunneling conductance, except that the DOS N((o) 
is replaced by an "effective tunneling DOS" Nt((o) = 
V (0)Re [\co\/ y / 0) 2 -A 2 {co)] . N (0) is the normal-state DOS 
at zero energy — Nq((o) = Nq(Q) is assumed — and A(ft)) is 
the gap function. The latter must be understood as a Fermi- 
surface average of weakly momentum-dependent quantities, 
A(o) = (4>(fc,fi))/Z(Jt,0)))Fs with 4> and Z the Eliashberg 
pairing and renormalization functions. In the notation of 
Eq. (2), they read Z(k, G)) = l-[±u (k, (O) + E22O, ca)}/{2(0) 
and, for an s-wave gap Ao, 4>(&, ft>) = Aq + T.n{k, ft)). In 
a if-wave superconductor, the Fermi-surface average of the 
gap Aj vanishes, and so does the average of the off-diagonal 
self-energy since £12(^,0) « A^. The effective tunneling 
DOS concept is logically generalized 19 - 20 by writing Nj{(o) = 

with 



Vo(0)Re(|ft)|/yft^ 



<Kft>) = 



[A^(to)] 2 ) FS 

-±n{k,(o)/A k 



1 



1 - pii(*,fl»)+Z22(t,fl»)]/(2fl») / FS 



(6) 



This form of Nj(co) is an even function of ft), and cannot fit 
the particle-hole asymmetric spectra in HTS. Therefore, a fur- 
ther generalization of the effective tunneling DOS has been 
necessary, namely, 



Nt(co) = Vo(ft))Re 1 



v/ft. 2 -[A^(ft.)] 2 



(7) 



FS 



which suggests that the "true" superconducting DOS can be 
obtained by dividing the tunneling spectrum in the supercon- 
ducting state by the spectrum in the normal state. 20 ' 45 ^ 7 

Equation (7) is very convenient, but lacks a formal justi- 
fication. Our model offers the opportunity to investigate the 
usefulness of Eq. (7), by comparing numerically the actual 
tunneling DOS N(co) of Eq. (5) with the effective tunneling 
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FIG. 6. (Left panels) Real part (solid lines) and imaginary part 
(dashed lines) of the pairing function defined in Eq. (6) for t z = 
(2D) and t z = 200 meV (3D). (Right panels) Comparison of the ef- 
fective tunneling DOS N T (a>) of Eq. (7) with the actual DOS N(co) 
of Eq. (5). The dashed lines show Nt((d)/Nq(cq). In all graphs, the 
dotted vertical lines mark the energy ±(£2 sr + Aq). 



DOS Nt((o). For the practical evaluation of Nt(co), we de- 
fine the Fermi-surface average as 



/FS 



E*Aq(*,0) (■■■) 
E*Ao(*,0) 



(8) 



with Ao(Jk,0) the zero-energy spectral function in the ab- 
sence of pairing: Aq(Ic,0) = (— l/?r)ImGi 0)^=0- With 
this definition, the average is performed on the renormalized 
Fermi surface, defined by £jj + Re £ 1 1 (k, 0) = 0, rather than the 
bare Fermi surface ^ = 0. Furthermore, each state gets cor- 
rectly weighted if the spectral weight is unevenly distributed 
along the Fermi surface. 

A comparison of N((o) and Nt(co) is displayed in Fig. 6, 
where the pairing function <p (co) is also shown. In two dimen- 
sions, the real part of (co) has a maximum at CO = Cl ST + A(>, 
where its imaginary part shows a rapid variation. This is 



analogous to the behavior reported in Ref. 13. The result- 
ing Nt(co) also shows a behavior similar to the one found 
in Ref. 13: Nt(co) is larger than the BCS density of states 
at energies smaller than il sr + Aq and drops below the BCS 
DOS at i2 sr + A . The actual DOS N(co), however, behaves 
differently: it is smaller than the BCS DOS between the co- 
herence peak and some energy above the dip minimum (see 
also Fig. 1 of Ref. 8). Thus, although the positions of the 
strong-coupling features are identical in Nt(co) and N(co), 
their shape is markedly different in 2D cZ-wave superconduc- 
tors. In 3D, the difference between Nt(co) and N(co) is less 
severe than in 2D, and both curves show very weak signatures, 
although those in N(co) are slightly stronger. Finally, the 
Nt(co) curves show structures which are absent in the N(co) 
curves. In 2D, a peak at co = — 16 meV = ^.o) appears due to 
the VHS in Nq(co); this peak is unphysical because in the ac- 
tual energy spectrum, the VHS is pushed to — [& n ^ + Aq] l l 2 . 
In 3D, Nt(co) has a structure near —100 meV, which also 
comes from the bare DOS Nq(co) as can be seen in Fig. lb. 
In the actual spectrum, this structure is suppressed due to the 
persistence of a large scattering rate at energies much higher 
than the threshold £2 S1 - (see Fig. 2). These problems illustrate 
the limitations of the simple product Ansatz Eq. (7) for ana- 
lyzing the tunneling spectrum of c/-wave superconductors. 
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APPENDIX: ANALYTICAL k z INTEGRATION 

If the Nambu self-energy has no k z dependence, the k z sum 
in Eq. (5) can be performed analytically. This is the case in our 
model defined in Eq. (2). Solving Dyson's equation G(k, co) = 
[Gq 1 (it, co) - l(k, co)}- 1 with Go given by Eq. (3), we find 



d i (*, co) ^ + , r _^_ tu ^ w j _ ^ + £ i2(jfe) ^2/ [co + iY + & _ t 22 (k 1 co)] ■ 

Since £ does not depend on k z (although it does depend on f-), the k z dependence only comes from ^ in Eq. (4) and we can 
make it explicit by rewriting 

G n (k,co) = 



zii -2f z cos(£ z c) -Zj 2 /[z22 + 2f-cos(A:-c)] 



A / Tj + A — 2f,cos(£ z c) 2 \ X) 77 -A -2f,cos(^c) 

In Eq. (10), the quantities zu> Z22, Zn, C> ^> an d ^ ^ ai l functions of ife|| and CO but not of k z . Explicitly, z\\ = 03 + iT— ^ — 
t n (k«,o), Z22 = co + rT+4 -£2 2 (fc||,Ci)), zn=\ +t n (k ll ,co), C = (zn +z 2 2)/2, A = J?-z\ 2 , and 77 = (zn -z 22 )/2. 
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The k z integration can then be performed by means of the and yields 
identity, 



l-Z-)D k {T]-X)}, (12) 



1 /*" dx 1 

b t (z) = — / = , (11) where JV\\ is the number of k\\ points in the 2D Brillouin zone. 
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